Structure of Flux Line Lattices with Weak Disorder at Large Length Scales 
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Dislocation-free decoration images containing up to 80,000 vortices have been obtained on high 
quality Bi2Sr2CaCu2 08+i: superconducting single crystals. The observed flux line lattices are in the 
random manifold regime with a roughening exponent of 0.44 for length scales up to 80-100 lattice 
constants. At larger length scales, the data exhibit nonequilibrium features that persist for different 
cooling rates and field histories. 

PACS numbers: 74.60.Ge, 75.10.Nr, 74.60.Ec, 74.72.Hs 



Recent studies of high temperature superconductors 
have shown richness in the phase diagram due to the pres- 
ence of weak quenched disorder |l| . Larkin first showed 
that arbitrarily weak disorder destroys the long range 
translational order of flux lines (FLs) in a lattice [g[. It 
was recently pointed out that the Larkin model, which is 
based on a small displacement expansion of the disorder 
potential, cannot be applied to length scales larger than 
the correlated volume of the impurity potential termed 
the Larkin regime P-p|. Beyond the Larkin regime, the 
behavior of FLs in the absence of dislocations has been 
considered using elastic models |^-|^. First, FLs start 
to behave collectively as an elastic manifold in a random 
potential with many metastable states (the random man- 
ifold regime) 0]. In this random manifold regime, the 
translational order decreases as a stretched exponential, 
whereas there is a more rapid exponential decay in the 
Larkin regime. At even larger length scales, when the 
displacement correlation of FLs become comparable to 
the lattice spacing, the random manifold regime transits 
to a quasiordered regime where the translational order 
decays as a power law |Q,g. 

Experimentally, neutron diffraction |6| and local Hall 
probe measurements 0] have shown the existence of an 
order-disorder phase transition with increased field, al- 
though the microscopic details of these phases are not 
clear. Theoretical progress describing FLs in the pres- 
ence of weak disorder has been made within elastic the- 
ory, whichproposes the absence of dislocations at equilib- 
|. To date, however, there has been no exper- 



imental work addressing the structure of dislocation-free 
FL lattices at large length scales. Previous magnetic dec- 
oration studies [ [lO|Jll| ] showed that the dislocation den- 
sity decreases and the translational order increases with 
increasing magnetic field. However, only relatively short- 
range translational order could be probed in the previous 
work due to the finite image size and relatively low ap- 
plied fields. 

In this paper, we report the first large length scale 
structural studies of FLs with measurements extending 
up to ~ 300 lattice constants and fields up to 120 G. Real- 
space images show dislocation free regions containing up 
to the order of 10^ FLs. A very low density of dislocations 



was also observed, although detailed analysis suggests 
that the dislocations are not equilibrium features. The 
translational correlation function and displacement cor- 
relator have been calculated from dislocation free data to 
examine quantitatively the decay of order. These results 
show a stretched exponential decay of the translational 
order indicating that FLs are in the random manifold 
regime. The experimentally determined roughening ex- 
ponent in the random manifold regime agrees well with 
theoretical predictions. 

High quality single crystals of Bi2Sr2CaCu208+2; 
(BSCCO) were grown as described elsewhere |12J. Typ- 
ically, crystals of ~ \mm x 1mm x 20/im size were 
mounted on a copper cold-finger and decorated with ther- 
mally evaporated iron clusters at 4 K. The samples were 
cooled down to 4 K using different thermal cycles to test 
nonequilibrium effects and to achieve as close an equi- 
librium configuration of FLs as possible within the ex- 
perimental time scale. The FL structure was imaged 
after decoration using a scanning electron microscope 
equipped with a 4096 x 4096 pixel, 8-bit gray-scale im- 
age acquisition system. Nonlinearity in the system was 
eliminated using grating standards. This high-resolution 
system enabled us to acquire images containing nearly 
10^ FLs, while maintaining a similar resolution (~ 14 
pixels between vortices) to previous studies of lO'^ FLs. 
In addition, an iterative Voronoi construction |13| was 
used to reduce the positioning inaccuracy to 3 % of a 
lattice constant. 

Samples were decorated in fields of 70, 80 and 120 
G parallel to the c axis of BSCCO single crystals. In 
contrast to the previous decoration experiments at lower 
fields [ [lO|jll]] , we find that the dislocations are rare at 
these fields. The density of dislocations was 1.7 x 10~^, 
1.4 X lO"-"^ and 3.1 x IQ-^ for 70, 80 and 120 G, respec- 
tively, where the total number of vortices is ^ 240, 000 for 
each field. It is thus trivial to find many large 100 x 100 
jiTT? dislocation free regions in the decorated samples. 
The size of the largest dislocation free image, which was 
obtained in a field of 70 G, is 152 x 152 ^itc? with 78,363 
vortices pj]. Although a small number of dislocations 
are detected in our FL images, this does not imply that 
they are energetically favorable at equilibrium. On the 



contrary, we believe that the large dislocation-free areas 
observed in the images provide a lower bound for the 
length scale of equilibrium dislocation loops. We discuss 
this point below after presenting a quantitative analysis 
of the translational order. 

To study quantitatively the FL lattice order, we pro- 
ceed as follows. First, a perfect lattice is constructed and 
registered to the FL positions obtained from an experi- 
mental image. The initial lattice vectors used to con- 
struct the perfect lattice were obtained from the Fourier 
transform of the vortex positions. When an image con- 
tains a dislocation, the continuum approximation is used 
to construct the perfect lattice with the dislocation [g5| . 
Next we minimized the root mean square displacement 
between the underlying perfect lattice and the real FL 
lattice by varying the position and orientation of the two 
lattice vectors of the perfect lattice. The displacement 
vector u(r) associated with each of the vortices posi- 
tioned at r relative to the perfect lattice was then com- 
puted. Fig |] displays a color-representation of the dis- 
placement field for a typical dislocation-free image and 
an image containing three dislocations. In Fig |l|(a), the 
average displacement is 0.22 oq, where ag is the lattice 
constant. Qualitatively, the map consists of several in- 
termixed domain-like structures, within which the dis- 
placement fields are correlated. These uniformly dis- 
persed domain-like structures of the displacement field 
produce sharp Bragg peaks in Fourier space (see Fig 0(b) 
later). We also believe that u(r) provides a quick indi- 
cation of nonequilibrium effects. For example. Fig 0(b) 
exhibits large domains of correlated displacements that 
are sheared relative to each other; that is, the blue-green- 
blue coded domains. We believe that this larger scale dis- 
tortion is a manifestation of a nonequilibrium structure 
that may arise from quenched dynamics of FLs during 
our field-cooling process (see below). 

To compare our data directly with theoretical predic- 
tions [g|, we have calculated the displacement correla- 
tor, B{r), and translational correlation function, Cci{r). 
Bir) and Cg(?') are defined as ([u(r) — u(0)]^)/2 and 
^giG[u(r)-u(o)]^^ j-ggpgctively, where ( ) is the average over 
thermal fluctuations and quenched disorder, and G is 
one of the reciprocal lattice vectors. Theoretically ^, 
we expect B{r) will show three distinct behaviors as r in- 
creases: B{r) '^ r in the Larkin regime, where B{r) is less 
than the square of ^, the in-plane coherence length. As r 
increases further, FLs are in the random manifold regime 
where ^^ < B{r) < Oq. In this regime B{r) r^ r^" with 
the roughening exponent 2i^ (< 1). Finally, at the largest 
length scales (the quasiordered regime) where a§ < B[r), 
B{r) ^ Inr. Since the in-plane ^ of BSCCO is only 
~ 20 A, the Larkin regime is irrelevant in our exper- 
iment (i.e., flQ » S,)- Fig 0(a) shows the behavior of 
B(r) calculated from the data in Fig 0(a). For r < 80ao, 
B{r) can be fit well with a power law, B{r) ~ r^", with 
21^ = 0.44. Thus our experiment is probing the ran- 



dom manifold regime at least up to this scale. Indeed, 
B{r) grows only up to 0.05 Og at r = SOoq, well below 
the expected crossover to the quasiordered regime, i.e. 



Bir) 



A naive extrapolation to B{r) — ag sug- 



gests the crossover at r ~ 10,000ao (~ 4 mm), which is 
far beyond our experimental limit. Samples with such 
a large clean area, and direct imaging of ~ 10® vortices 
would be required to observe the logarithmic roughening 
of FLs. The roughening exponent 2u is found to be inde- 
pendent of the field (70 - 120 G) and consistent with the 
estimate 2^ = 2/5 obtained by Feigelman et al. using a 
scaling argument S. As shown in Fig 0(b), Ccir) and 
g~G B(r)/2 overlap with each other for r < L*, where the 
measured L* is ~ SOag. These results support the Gaus- 
sian approximation, Ccir) « e~'^ B{r)/2^ which has been 
simply assumed for the equilibrium FLs lattice B within 
this length scale. For r > L*, however, B(r) deviates 
strongly from expected behavior; that is, B{r) saturates 
and even decreases as r increases. In addition, the Gaus- 
sian approximation breaks down for r > L* as evidenced 
by the difference between Ccir) and e~^ B(r)/2^ -y^g 
believe that this behavior can be attributed to nonequi- 
librium FL structures at the larger length scales of our 
experiment. 

To examine this point further, we decompose B{r) 
into its longitudinal [i?'^(r)] and transverse [i3"^(r)] parts: 
B{r) = {B^{r) + B'^(r))/2, where 



(1) 



B^(r) = ( (u(r)-u(O)).- ). 



It is worth noting that in the random manifold regime, 
the ratio of B^{r) and B^{r) is predicted to be 2z^-f 1 [|1, 
and thus an independent estimate of the roughening ex- 
ponent. The average value of this ratio measured from 
our data (inset to Fig |^(a)) is 1.40, which is consistent 
with the value of 2v obtained from B[r). As shown in 
Figy(a), both i?^ and B^{r) are described well with the 
power law behavior up to r '^ L*. Beyond this range, 
however, the transverse displacement B^{r) first devi- 
ates from power law causing deviations in B{r). Thus, 
we infer that shear motion of FL lattice should be re- 
sponsible for the abnormal behavior of B{r). Since the 
shear modulus of FL lattice is much smaller in magni- 
tude than the compressional modulus ||l|, B^{r) is al- 
ways larger than B^ (r) , and the shear motion dominates 
the relaxation of the FL lattice during the field cooling 
process. As temperature decreases, the long wavelength 
component of shear motion is frozen out. We believe 
that the domain-like structures seen in Fig ^ are a snap 
shot of these frozen long wavelength shear motions. Note 
that the characteristic length scale of these domain like 
structures in Fig 1 is again ~ L*, which explains the 
deviations in B{r) for r > L* . Therefore, L* is the equi- 
librium length scale within which FLs can relax to the 
local equilibrium during our experimental time scale. 



This issue can also be addressed through Fourier space 
analysis. Fig ||(b) displays a blow-up of one Bragg peak. 
Several small satellite peaks appear around the relatively 
sharp main peak; these satellite peaks indicate a large- 
scale modulation of the FL lattice. If the FLs were in 
equilibrium, only one main peak should be expected. The 
corresponding real space distance between the main and 
satellite peaks is, again, ~ L*. Hence these satellite peaks 
provide another evidence of the frozen-in dynamics be- 
yond the equilibrium length scale L*. In addition, we 
have prepared FL lattices in different ways to address 
the nonequilibrium structures. For example, we cooled 
the samples in the absence of a field to 65 K, applied a 
field 70 G, and then cooled slowly (0.1 K/min) to 4 K. 
Significantly, we find a similar density of dislocations and 
FL structure compared to the rapid (10 K/s) field-cooled 
samples. Since 65 K is far below the melting tempera- 
ture |l^, this observation suggests that the nonequilib- 
rium structures originate from the frozcn-in dynamics far 
below the melting temperature. Although we can probe 
FLs up to a length scale of ~ 300 uq, there is a much 
smaller length scale L* that prohibits direct application 
of the theory derived for an equilibrium FLs. Further 
studies should address this important issue. 

Finally, we consider the origin of dislocations observed 
in our experiments, since nonequilibrium vs. equilibrium 
nature of dislocations is critical to the existence of the 
Bragg glass phase. We believe that our data, which ex- 
hibit the small numbers of dislocations, in fact, favors 
nonequilibrium nature of dislocation in the FL lattice we 
probed by following reasons. First, it is found that most 
dislocations are pinned in between domain boundaries 
(see Fig WPo) for example). If there were a dislocation 
within the domain-like structures where FLs are locally 
in equilibrium, the dislocation should be an equilibrium 

— 1/2 

feature. Second, L* <^ Ld = n^ ' '^ 250ao, where 
Ld and n^ are the average distance between dislocations 
and the density of dislocations, respectively. If disloca- 
tions were energetically favorable in an equilibrium FL 
lattice, large dislocation loops should proliferate beyond 
the equilibrium length scale L* . In addition, if some dis- 
locations drift within domains, and are pinned at domain 
boundaries, we should have Ld ^ L*. Therefore, our ex- 
periment {L* 4C Ld) suggests that dislocations are not 
equilibrium features in the FL lattice. Together, our data 
provide a lower bound for the length scale of equilibrium 
dislocation loop in the FL lattice. 

In summary, we have obtained large scale dislocation- 
free images of the FL lattice in high quality BSCCO su- 
perconductors. Quantitative analyses of the translational 
order indicate that the system is in equilibrium for length 
scales up to ~ SOoq, and that FLs are in the random man- 
ifold regime with a roughening exponent 2i^ = 0.44. We 
suggest that the very small density of dislocations ob- 
served in our data is an out-of-equilibrium feature due to 



the short time scales involved in our field-cooled experi- 
ments. 
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FIG. 1. (color) Spatial map of vortex displacements u(r) 
from the perfect lattice positions. Darker (brighter) regions 
in the map correspond to smaller (larger) displacements. Dif- 
ferent colors correspond to vortex displacements in different 
directions, as shown in the inserted color wheel. The two 
solid lines, inner and outer circles in the color wheel corre- 
spond to two basis vectors of the lattice, displacements of 0.5 
ao and ao respectively. Samples were decorated at 70 G. (a) 
dislocation-free image containing 37003 vortices. The edge of 
the image correspond to 106 /im. The lower inset shows a 
part of both real FL image and a perfect lattice (yellow) with 
displacement vectors (red), (b) larger scale image containing 
three dislocations (highlighted by red dots and circles) and 
large scale shearing. The image contains 78385 vortices in 
the 160 X 160 nm^ area. 



FIG. 2. (a) Mean-square relative displacement correlator 
B{r) (open circles) as a function of distance r calculated 
from the image shown in Fig 1(a). The solid line is a power 
law fit : B{r) ~ r'''*''.(b) Translational correlation function 
Cg('") (dots) calculated from the same image. The open 
circles are a comparison with the Gaussian approximation : 



FIG. 3. (a) Transverse (open circle) and longitudinal 
(solid circle) displacement correlators, i.e. B^ {r) and B^{r), 
as a function of distance calculated from Fig. Ma). The insert 
shows the ration of the two quantities, (b) The Fourier space 
image showing all six first order Bragg peaks, (c) Detailed of 
one of the Fourier peaks calculated from the same image (in- 
verted gray scale). Dark arrows highlight two satellite peaks. 



This figure "Fig 1. GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/99 1 1 84v2 



This figure "Fig2.GIF" is available in "GIF" format from: 
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This figure "Fig3.GIF" is available in "GIF" format from: 
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